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1. Introduction 

The phase diagram of QCD in the region of vanishing chemical potential is actively explored in 
heavy ion collision experiments at RHIC and LHC. Lattice calculations (see [1] for a recent review) 
predict a crossover transition to a larger number of degrees of freedom, typical of a deconfined 
medium (the Quark-Gluon Plasma), for temperatures T ranging from 150 to 200 MeV. On the 
experimental side, the characterization of the properties of the medium relies either on its bulk 
properties, which find an effective description in hydrodynamics, or in hard probes, i.e. energetic 
particles not in equilibrium with the medium. 

Heavy quarkonium has been one of the most actively investigated hard probes for the past 26 
years. In 1986 Matsui and Satz [2] hypothesized that colour screening in a deconfined medium 
would have dissociated the J/i\f, resulting in a suppressed yield in the easily accessible dilepton 
decay channel and yielding a striking QGP signature. Such a suppressed yield has been indeed 
observed in heavy ion collision experiments at SPS, RHIC and LHC (see [3] for a review), the 
current understanding being that it cannot be explained by cold nuclear matter effects alone, i.e. 
those caused by a confined nuclear medium. Furthermore the LHC results have opened up the 
frontier of the cleaner bottomonium probe with the availability of quality data on the T resonances 
(see [4] for the latest CMS results). 

On the theory side, a great deal of the studies of the in-medium dynamics of the QQ bound 
states has been carried out with phenomenological potential models, fust introduced in [5], where 
all medium effects are encoded in a T-dependent potential plugged in a Schrodinger equation. 
We refer to [6 - 8] for recent reviews. The derivation of such models from QCD was however 
not established. Moreover, lattice calculations of free energies and other quantities [9, 10] obtained 
from correlation functions of Polyakov loops are often taken as input for the T-dependent potential. 
Although these quantities have been thought to be related to the colour-singlet and colour-octet 
heavy quark potentials at finite temperature [11, 10], a precise relation was still lacking in the 
literature, as pointed out in [12]. 

In the past few years, a considerable effort has been devoted to the derivation of a proper 
potential from QCD. In [13] (see also [14]) it was shown that the perturbative real-time potential 
obtained by analytical continuation of a rectangular Wilson loop shows a screened real part and a 
large imaginary part (for T S> 1/r ~ nip, mp being the Debye mass, the imaginary part is much 
larger than the real part) which is due to the scattering of the virtual potential gluons with the 
light particles (quarks and gluons) of the medium, in what is called Landau damping. Further 
phenomenological studies have shown the numerical importance of this imaginary part [15-17] 
when solving a Schrodinger equation. A generalization to anisotropic plasmas has been performed 
in [18-20] and used in [21, 22] (coupled with a hydrodynamical evolution) to obtain predictions 
on the suppression rates. 

The first-principle derivation of potentials from QCD was further improved with the develop- 
ment of an Effective Field Theory (EFT) framework [23-31], which will be at the center of this 
contribution. This approach, started in [23] (see also [24] for the abelian analog), relies on the 
well-established framework of T = Non-Relativistic (NR) EFTs of QCD and extends it to finite 
temperatures, giving to the potential a modern, rigorous definition as a matching coefficient of the 
EFT that arises when the scale mv is integrated out, where m is the heavy quark (Q) mass and v <C 1 
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is the relative velocity of the constituents of the NR bound state. At zero temperature this EFT is 
called Potential Non-Relativistic QCD (pNRQCD) [32, 33]. 

This contribution is structured as follows: in Sec. 2 we will illustrate the basic principles of 
the EFT framework. In sections 3 and 4 we will apply those principles to the real-time dynamics 
of bound states in two interesting temperature regimes, T > mv and T < mv respectively. In Sec. 5 
we will compare the thermal widths obtained within the EFT framework with other approaches in 
the literature and show how dissociation cross sections can be extracted from the EFT calculations. 
In Sec. 6 we will deal with Euclidean quantities and their related thermodynamic al free energies, 
showing how they relate to the QQ potentials. In Sec. 7 we will briefly review some interesting 
recent developments and finally in Sec. 8 we will draw our conclusions. 

2. EFTs 

In a scenario where a physical system is described by two scales A»m,a low-energy EFT 
describing the physics at the scale m is constructed by integrating out A, i.e. integrating out modes 
of energy /momentum above a given cut-off jU, A S> fL 3> m, and imposing that the EFT and the 
theory it is derived from are equivalent below the cut-off. If the system is characterized by more 
than two scales and if these scales are hierarchically ordered, the procedure can be iterated, creating 
a tower or hierarchy of EFTs. If two scales are not hierarchically separated, they should then be 
integrated out in the same step. 

This line of action is at the base of the above-mentioned T = NR EFT framework. As 
we remarked, the expansion parameter there is v <C 1, giving rise to the hierarchy of scales m S> 
mv S> mv 2 , where m is the heavy quark mass, mv is the typical relative momentum or inverse 
radius and mv 2 the kinetic/binding energy. The framework is constructed by integrating out first 
the mass scale, yielding Non-Relativistic QCD (NRQCD) [34, 35]. Its Lagrangian is organized 
as an expansion in the inverse mass. The second step is the integration of the scale mv, which 
leads to pNRQCD [32, 33]. The relative hierarchical position of Aqcd and mv establishes whether 
this integration is to be performed perturbatively or non-perturbatively. The former case leads 
to weakly-coupled pNRQCD, which we illustrate briefly because the EFTs developed at finite 
temperature, which will be introduced later, will share its basic form. 

The degrees of freedom in the QQ sector are conveniently written as a colour-singlet and a 
colour-octet QQ bilinear field, which can interact through ultrasoft (E,p~ mv 2 gluons). This cutoff 
is effectively expressed through a multipole expansion. In detail, the Lagrangian reads [32, 33] 

1 " f r r 

^pNRQCD = - t V" V + £ ?i iP U + / d 3 r Tr { S t [id - h s ] S + O 1 " [iD - h ] O 
4 i=i J L 

+V A (O t r^ES + H.c.) + ^O t {r- 5 E,0} + ...}. (2.1) 

The fields S = S\ c /\fNc and O = 0"T" /\/Tp are the QQ colour-singlet and colour-octet fields 
respectively, qi are. the light quarks, in nj flavours, Tp = 1/2, E is the chromoelectric field, iD$0 = 
idoO — gA()0 + OgAo and H.c. stands for Hermitian conjugate. The trace is over colour and spin 
indices. The dots in the last line stand for higher-order terms in r and 1/m. 
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The dependence on the scales m and mv ~ ma s is encoded in the Wilson coefficients; Va and 
Vg are at leading order Va = Vg = 1, whereas the singlet and octet Hamiltonians have the form 
(P = -f'V r ) 

h s = C F h..., K = h— V-... (2.2) 

m r m 2N C r 

The dots in these equations stand for higher orders in the 1/m and a s expansions, such as spin- 
and momentum-dependent terms or radiative corrections to the Coulomb potentials. N c = 3 is the 
number of colours and Cf = 4/3 in QCD. In the power-counting of the EFT the explicitly-shown 
terms are of the same size and give rise to a Coulombic bound state. At the zeroth order in the mul- 
tipole expansion the equation of motion for the singlet field arising from Eq. (2.1) is then a simple 
Schrodinger equation, resulting in a spectrum made by the (QCD) Bohr levels E n = — 4ma 2 /(9« 2 ), 
whereas the octet potential is repulsive and does not support bound states but a continuum of scat- 
tering states. Higher orders in the multipole expansion, such as the chromoelectric dipole couplings 
on the second line of Eq. (2.1), are responsibile for retardation effects such as those yielding the 
Lamb shift in hydrogen spectroscopy. 

At finite temperature one also encounters the scales that characterize the thermal medium: its 
temperature T, the electric (Debye) screening mass nip ~ gT and a magnetic mass of order g 2 T } 
In the weak-coupling scenario that has so far been investigated in the EFT approach, these scales 
develop a hierarchy. In order to construct a finite temperature EFT framework analogous to the 
one just illustrated, one needs then to put these scales together with the NR ones in a global hierar- 
chy. Many such hierarchies are possible and phenomenologically sensible. The latter requirement 
translates into r«m. 

Once a hierarchy is fixed, one proceeds to integrate out scales along the procedure sketched 
above. Since m 3> T, the first step is always the integration of the mass, leading again to NRQCD. 
Since an EFT and its matching coefficients are independent of the hierarchy of scales below the 
cutoff, its Lagrangian is identical to the zero-temperature case. 

After this first step, the subsequent ones differ in the different cases, which in the following 
we will group into two regions: the one where T 3> mv and the one where T <C mv, which can be 
of relevance for the ground state of bottomonium. The case T ~ mv has been dealt with in QED in 
[27]; its QCD investigation is underway in [36]. When integrating out the thermal scales, existing 
finite-T EFTs of QCD such as the Hard Thermal Loop (HTL) effective theory [37] are employed. 

In all cases, once the scale mv is integrated out, the resulting Lagrangian is similar to the 
one of pNRQCD, with colour-singlet and colour-octet fields and a Schrodinger equation picture 
appearing as the zeroth-order equation of motion. The potential is then rigorously defined as a 
matching coefficient; all scales above mv will contribute to it. 

The different possibilities for the hierarchies and the EFTs that arise out of them are pictured in 
Fig. 1. We conclude this Section by noting that the NR EFT framework inherits all the strong points 
of EFTs. In particular, it has a well-defined power counting, which allows for easy estimations of 
the size of possible contributions to interesting observables, and it is systematically improvable, 

'We do not distinguish between T and TtT or multiples thereof in the text. For what concerns the magnetic mass, 
we never reach an accuracy where its contribution becomes relevant, so we need not worry about its non-perturbative 
treatment. 
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Figure 1: Possible hierarchies and the EFTs they generate. Figure taken from [38]. 

i.e. new operators and matching coefficient can be added to the Lagrangian and computed in order 
to reach the desired accuracy in the calculation of the above-mentioned observables. 



3. The m>r»mv region 

In this region the temperature scale is to be integrated first, obtaining an EFT of NRQCD called 
NRQCDhtl [39, 40], where the gauge and light-quark sector is described by the HTL Lagrangian. 
The second step depends on the existence of a hierarchical separation between mjj and mv. If 
nip ~ mv the two scales are to be integrated out at the same time. If, on the other hand, mv 3> mp, 
one needs two separate steps. For reasons that will become clearer later, the case mv <C mo is not 
relevant. 

In the first case one obtains a pNRQCD-like EFT named pNRQCDhtl, where the potential, 
obtained by matching from NRQCDhtl in the real-time formalism of Thermal Field Theory, re- 
produces the earlier result of [13], thus bringing it in a consistent EFT picture [23]. This derivation 
furthermore shows the computational advantage of the real-time formalism over the imaginary-time 
one for this class of Minkowskian observables. The colour-singlet potential reads 



V s (r)\ 



-C F — e~ nDr - C F a s m D + i2C F a s T 
r 



dt 



sin(«i£)rf) 
mort 



1 



(3.1) 



g T (1 +«//6). We remark that the real part scales like a s mo (recall that mv 



where m 2 D 

1/r ~ mo). The imaginary part is Bose-enhanced and scales like a s T; it is thus much larger, as we 
mentioned before. Hence, one can expect bound states to be dissociated in this regime. 

If one moves to the second possibility mv 3> mo, the integration procedure has to go through an 
intermediate step. The real part is dominated by the T = Coulomb term, whereas the imaginary 
part at the scale mv is IR divergent and is regulated in dimensional regularization. Once the UV- 
divergent contribution of the smaller scale nip is added up, the following finite result arises [23] 



VJr) 



-C F — - -^-a s rm 2 D + i-^-a^Trnj) ( \n{rm D ) £ + 2y E 



(3.2) 
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where Je is the Euler-Mascheroni constant. As pointed out in [41] (see also [24] for QED), the real 

2/3 9 

and imaginary parts of this potential become parametrically of the same size when T ~ ma s . 
This then constitutes a parametric estimate of the dissociation temperature, where the dissociation 
criterion popularized by potential models, i.e. disappearance of the binding energy, is replaced by 
the equality of binding energy and thermal width. A numerical estimate has been performed in [27] 
for the T(1S) and is reported in Table 1. 



m c (MeV) 


T d (MeV) 


oo 


480 


5000 


480 


2500 


460 


1200 


440 





420 



Table 1: Dissociation temperature for the T(15) for different values of the charm mass. Table and results 
are taken from [27]. 



4. The mv>T> mv 2 region 

In this region the thermal medium acts as a perturbation to a Coulombic bound state, but it still 
modifies the potential. Furthermore, as argued in [42], hierarchies where mv > T > mv 2 could be 
of relevance for the ground states of bottomonium, where mv ~ ma s ~ 1.5 GeV and T or even kT 
are smaller in current experiments. 

In this region one starts by integrating out mv from NRQCD, obtaining standard pNRQCD. 
The next step, the integration of the temperature, leads to a different version of pNRQCDhtl than 
the one considered above. Its Lagrangian has been written and some of its matching coefficients 
have been computed in [25, 29]. 

Further steps depend on the hierarchy (if any) between mv 2 and mo- In [25] the case mv 2 S> mo 
has been explored in detail, leading to a calculation of the thermal modifications to the spectrum 
and of the thermal width up to order ma s 5 in the power-counting of the EFT. Before we briefly 
illustrate these results, we remark that other hierarchical possibilities have been considered in the 
static case in [23] and in QED in [27]. 

In order to obtain the spectrum and width to order wia s 5 , as well as to compute the matching 
coefficients of pNRQCDhtl. one has to compute loop diagrams in pNRQCD with loop momentum 
of order T or in pNRQCDhtl with loop momenta of order mv 2 or mj). These loop diagrams arise 
from the dipole interaction vertices in Eq. (2.1); the relevant diagrams for the desired accuracy are 
shown in Fig. 2. As in the previous section, we encounter intermediate UV or IR divergences which 
we regulate in dimensional regularization. The potential, which in the EFT framework is not an 
observable, shows IR divergences. The spectrum and the width, being instead observables, are of 
course finite, the divergences cancelling in the sum of the contributions from different scales. 

2 The Coulombic real part implies a Coulombic power counting, with mv ~ l/r ~ mcc s . Logarithms of the coupling 
are not considered in this estimate 
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Figure 2: Loop diagrams in pNRQCD and pNRQCDhtl contributing to our calculation. Single lines are 
colour-singlet QQ states, double lines are colour octets, curly lines are gluons, vertices are chromoelectric 
dipoles and the blob is the gluon self-energy. The imaginary part of the first diagram yields the singlet-to- 
octet decay mechanism, whereas the second one gives the Landau damping contribution to the width. 



The shift to the binding energy induced by the thermal medium, which translates into a mass 
shift for the bound state, reads to order ma s 5 

34 



8E 



(thermal) | 
15 



Imv3>rs>mv 2 



27 



K a 2 T 2 a Q 



+- 



7225 / 47T 2 r 2 
324 l0§ " 



128 



+3a^a s r { -C(3)<%r 



-27b 



; m D ( 



(4.1) 



where ao = 3/(2raa s ) is the Bohr radius and Lis is the QCD equivalent of the Bethe log in the 
hydrogen spectrum. We refer to [43, 44] for its definition and evaluation, which yields L\$ = 
-81.5379. 

The thermal width reads instead 



-.(thermal) i 
IS 



1156 



.2 



;i 



a:T + 



7225 



162 



Eia: 



3 

32 



m 2 D 



E 2 

ln^+2 7£ 



3-log4-2 



C(2) 



+ 87rln2 a,T' 



+ — a s Tm 2 D alh s 



(4.2) 



where lis = —0.49673 is another Bethe logarithm. We refer to [25] for details. The terms on 
the first line of this equation arise from the imaginary part of the diagram on the left of Fig. 2 
and correspond to the absorption of a thermal gluon by the bound state and its subsequent decay 
into a colour-octet state, a process we call singlet-to-octet decay. The terms on the second and 
third lines, on the other hand, arise from the imaginary part of the second diagram and are instead 
due to scattering of the virtual gluon with the light constituents, i.e. to Landau damping. In the 
power counting of pNRQCDhtl singlet-to-octet decay is dominant over Landau damping as long 
as mv 2 > mo, as in this case. We will return to these processes in the next section. 

We conclude this section by reporting on the findings of [45, 46] (see also [47]). There, the 
authors employed the lattice formulation of NRQCD to access bb bound states on the lattice at 
finite temperature. In [46] they used Bayesian MEM techniques to reconstruct the spectral function 
and thus extract 8E and T as a function of the temperature. These results are shown in Fig. 3. 
The authors then proceeded to compare their results with the leading term of the EFT results, 
represented by the first term on the r.h.s. of Eqs. (4. 1) and (4.2). They found them to be compatible 
for nib = 5 GeV and a s w 0.4. 
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77T 
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Figure 3: Plots of T and 5£ in the vector (T) channel from the lattice calculation of [46]. Figures taken from 
that reference. The dashed line on the right is the leading term of our EFT calculation [25], |^7TOC s 2 T 2 ao, for 
m b = 5 GeV and a s = 0.4. 



5. The thermal width in the EFT and in the literature 

In the previous sections we have identified two processes that are responsible for the appear- 
ance of the thermal width: these are the singlet-to-octet decay (g + x ¥ — > {QQ)%, where is the 
bound state and {QQ)% the QQ pair in the colour-octet state) and Landau damping + — > 
P + (68)8> where p = g,q,q is a light parton). Analogous processes had been previously con- 
sidered in the literature, going under the names of gluo-dissociation [48, 49] and dissociation by 
inelastic parton scattering [50, 51] (see [52-54] for recent applications) respectively. In these ap- 
proaches a T = cross section for the process at hand, sometimes complemented by momentum- 
independent thermal masses, was convoluted over a thermal distribution for the incoming parton, 
i.e. 

r= / ^f p (q)a(q)v reU (5.1) 

where f p is either the Fermi-Dirac or Bose-Einstein distribution, q is the incoming momentum, 
with kinematic minimum for dissociation q m i n , and v re i is the relative velocity. Cross sections 
employed in the literature have been the Bhanot-Peskin one [55] for gluo-dissociation and [56] or 
more recently [57, 58] for dissociation by inelastic parton scattering. 

In [30, 36] (see also [59]) we have undertaken a comparison between the EFT results and these 
phenomenological approaches. For what concerns gluo-dissociation, in [30] we have shown that 
our results (such as the first line of Eq. (4.2)) correspond to the formulation of Eq. (5.1) and that the 
Bhanot-Peskin cross section of [55] represents a large-A^. limit of the one that can be extracted from 
the EFT framework, which in turn correctly incoiporates the effect of the repulsive octet potential 
on the final state. An equivalent cross section has also been derived in [60] (see also [61, 62]). The 
EFT approach and its power counting clearly constrain the validity region of the gluo-dissociation 
approach to temperatures smaller than mv, so that the dipole approximation employed is valid, and 
larger than mv 2 , the threshold for dissociation. If the screening mass becomes of the same order 
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or larger than this latter scale, Hard Thermal Loop resummation becomes necessary and the cross 
section is no longer temperature-independent. 

For what concerns dissociation by inelastic parton scattering, our analysis [36, 59] shows that, 
for processes with a light particle in the final state, a factor of 1 ±f p (q out ) should be added to 
Eq. (5.1), where the plus sign applies to the case of the Bose distribution (Bose enhancement) 
an the minus for the Fermi distribution (Pauli blocking). Our calculation differs from the one 
in [50, 51], where the cross section is identified with two times the one for heavy quark-parton 
scattering computed in [56], supplemented by momentum-independent thermal masses. In our 
case the Landau damping terms in the EFT approach encode the square of the diagrams for heavy 
quark-parton scattering, of those for heavy antiquark-parton scattering and the interference term, 
which depends on the properties of the bound state. We show that in the power counting of the 
weakly-coupled EFT this last term can never be neglected and hence the use of two times the 
simple heavy quark-parton scattering cross section is inappropriate at weak coupling. Our cross 
section furthermore does depend on the temperature in all regimes. We plot an example in Fig. 4 
for the process g + *P(lS) — > q + (QQ)%. The different curves are the cross section computed in 
different regimes, which show the validity of the employed approximations. 




Figure 4: The cross section for q + ^VilS) — > q+ {QQ)s [36, 59], normalized by the constant a cq = 
yinrifQs fl o/3' where «o = 3/(2ma s ) is the Bohr radius, with moat) = 0.1. The horizontal axis is the 
momentum of the incoming quark in units of ciq. The continuous black line is the cross section for T ~ mv, 
the dashed blue line is the result for mv T, mo mv 2 , obtained in the multipole expansion, and the dash- 
dotted red line is the result for T mv ^> mo- The last two curves describe the limiting cases for the first 
one, which smoothly interpolates between them. We recall that typical thermal constituents of the plasma 
have momenta of order T and hence T ~ mv translates into qao ~ 1, whereas the blue curve is to be trusted 
for qao <C 1 and the red one for qao >• 1 . 

6. Euclidean quantities 

As remarked in the introduction, potential models have often used thermodynamical free ener- 
gies measured on the lattice as input. Popular choices have been the gauge-invariant colour-average 
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free-energy, derived from the correlator of two Polyakov loops [9], and the colour-singlet free en- 
ergy [10], which is gauge dependent. 

In [26] the correlator of two Polyakov loops has been analyzed in detail. A short-distance 
perturbative calculation was carried out to NNLO. Euclidean pNRQCD was then employed to 
show how the correlator can be expressed in terms of gauge-invariant colour-singlet and colour- 
octet free energies, which have been computed at short distances. An important finding is that 
these free energies are not equal to the real-time potentials computed in the same temperature 
regime: not only are the free energies purely real, in contrast with the complex potentials, but the 
real parts of the potentials are themselves different from the free energies. The discrepancy has 
been traced back to the different boundary conditions in the two cases: for the potentials one takes 
a large Minkowskian time to obtain the pole of the bound-state propagator, whereas for the free 
energy one has the Euclidean time spanning the entire compactified axis. 

The effect of periodic boundary conditions on another gauge-invariant quantity has been an- 
alyzed in detail in [31] (see also [63]). There, the renormalization of cyclic Wilson loop, i.e. a 
rectangular Wilson loop spanning the entire Euclidean time axis, has been studied. This operator 
had been found to be UV divergent after charge renormalization in [64]. The analysis of [3 1] shows 
how this extra divergence is due to the periodic boundary conditions, and how its renormalization 
involves a mixing with the Polyakov loop correlator. In particular, the difference W c — P c , W c being 
the cyclic loop and P c the correlator of Polyakov loops, is multiplicatively renormalizable and could 
then be another interesting quantity to be measured on the lattice and compared with perturbative 
result, in order to test the applicability of perturbation theory to <2<2-related quantities. 

7. Non-static and non-perturbative aspects 

The calculations reported in sections 3 and 4 are for a bound state at rest in the medium 
frame. In [28] (see also [59]) the NR EFT approach was extended to finite relative velocities 
between the medium and the (abelian) bound state. All the available velocity range, from v > 
to v — > 1, was explored. For muonic hydrogen, which is the abelian bound state more closely 
resembling quarkonium in a medium, it was found that, in the ultrarelativistic limit and for T > 
mv, the imaginary part due to Landau damping loses its role of main dissociation agent in favour 
of colour screening, whereas for smaller velocities, < v < 0.9, the bound state is more easily 
dissociated for increasing velocities. 

The effect of a finite relative velocity was also explored in lattice NRQCD in [65] (see also [47]) 
for bottomonium 5-wave states. A momentum dependence of the Euclidean correlators and the 
reconstructed spectral functions was observed and compared to the EFT results. 

A second topic of interest and focus of recent research is the non-perturbative determination of 
the potential. In [66, 67] the spectral function of the Wilson loop has been determined by applying 
the Maximum Entropy Method to lattice data and a complex potential has then been extracted 
from it. A different method has been used in [68, 69], where a real quantity has been fit from 
the exponential decay of correlators of temporal Wilson lines of varying time extent at a fixed 
spatial distance. The strong-coupling version of pRNQCD [33, 70] could be generalized to finite 
temperatures along the way the weakly-coupled one has; this would help shed more light on how 
the proper in-medium potential is to be non-perturbatively determined. 



10 



Review of the EFT treatment of quarkonium at finite temperature 



Jacopo Ghiglieri 



8. Summary and conclusions 

In this contribution we have briefly illustrated the development and a few chosen results of 
the NR EFT framework for in-medium quarkonium bound states. From the more theoretical point 
of view, we have seen how with EFTs one can give a rigorous QCD definition and derivation of 
the potential, bridging the gap with potentials models which appear as leading-order picture in 
pNRQCD and its finite-temperature analogues. Furthermore, the power counting EFTs naturally 
possess allow one to systematically take into account corrections, be it relativistic 1 /m n effects 
(n> 1), radiative corrections, etc. and include all medium effects which are relevant for the desired 
accuracy. 

For what concerns phenomenology, we have shown how previous results, such as the complex 
potential of [13], are accommodated in the EFT and how new results are obtained [23], such as 
the estimation of the T(15) dissociation temperature [27] given in Sec. 3 or the calculation of the 
corrections to the potential, the spectrum and the width for a Coulombic bound state in the regime 
mv>7'> mv 2 [25], illustrated in Sec. 4, which could be of relevance for the new frontier of T 
phenomenology in LHC experiments. We have also shown how the thermal width that arises in 
the EFT calculations compares with those obtained in phenomenological approaches based on the 
convolution of T = cross sections. For what concerns gluo-dissociation, the factorization formula 
holds and the cross section [55] is a large-A^ limit of the EFT cross section [30]. For dissociation 
by inelastic partem scattering, instead, our factorization formula differs from the one employed in 
the literature and a direct comparison between the cross sections is not possible due to different 
approximations and validity regions [36]. In both cases the power counting allows to constrain the 
validity region of these approaches. 

The EFT framework allows also an analysis of the thermodynamical free energies extracted 
from correlators of Polyakov loops and widely employed as input for potential models. It is found 
that the colour-singlet and colour-octet free energies that can be defined in the Euclidean EFT 
framework differ from the real part of the Minkowski-time potentials, the difference being traced 
back to the different boundary conditions in the two cases [26], as summarized in Sec. 6. There 
it is also shown how periodic boundary condition can induce UV divergences in these Euclidean 
operators, which then mix under renormalization [31]. 

Possible extensions of the framework include its generalization to the strong-coupling regime, 
which could help determining which Euclidean operators, amenable to lattice measurements, can 
be used to extract, possibly through analytical continuation, the real-time potentials governing the 
evolution of the bound state. 
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